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Abstract: We investigate a system described by two real scalar fields coupled with gravity 
in (4, 1) dimensions in warped spacetime involving one extra dimension. The results 
show that the parameter which controls the way the two scalar fields interact induces the 
appearence of thick brane which engenders internal structure, driving the energy density 
to localize inside the brane in a very specific way. 
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1. Introduction 

This work deals with scalar fields coupled to gravity in warped spacetime in (4,1) dimen- 
sions. We focus mainly on generation of thick branes, their stability and internal structure, 
and our investigations follow several recent works on the subject jl], [2], ||, [|, |5], ^, [?], ||, 



The brane world scenario that we investigate is described by scalar matter coupled 
to gravity. Thus, in flat spacetime the scalar matter model reduces to a system of scalar 
fields which supports interesting domain wall solutions. As one knows, domain walls are 
extended structures which may be seem as interfaces separating two distinct phases of the 
system. They appear in several branches of nonlinear science, in particular in high energy 
physics and in condensed matter |l7| , |i"8| . In condensed matter, for instance, domain 
walls may spring in several different ways, depending on the richness of the systems under 
consideration. An interesting situation occurs in ferromagnetic systems, as reported for 
instance in [17, 15], in which one finds walls of both Ising and Bloch type. Ising walls 



are simpler interfaces, having no internal structure. However, Bloch walls are interfaces 
that engender internal structure, and for this reason they may be seen as chiral interfaces. 



The chiral feature of Bloch walls was first explored in Ref. [19]; other recent investigations 
on Bloch walls in parametrically driven systems may be found in |2C], and in references 
therein. 

In high energy physics, domain walls usually appear in models described by real scalar 



fields [16]. In particular, in Ref. [pi] ] one introduces a simple model, described by two real 
scalar fields, which support both Ising and Bloch walls - see, e.g., Ref. [^, |23|, ^] for further 
details on the two-field model, and on Ising and Bloch walls. The approach introduced in 
Ref. [21] engenders an important issue, which concerns the Bogomol'nyi bound |25]. In 
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Ref. 21, 22, ^3|, 24 1 it was shown that this bound can be obtained in two-field models 
defined by the potential 



V(<f>,x) 



1 fdW\ 2 1 fdW x 2 
+ 



2 \d X 



where W = W((f>, x) is a smooth function of the two fields 
described by (^T) 

Wrfax) =<l>-\<t?-r<t>X 2 



(1.1) 



and x- The two-field model is 



(1.2) 



where r is a real parameter, and we are using dimensionless fields and spacetime variables. 
The richness of the two-field model was further explored in other contexts in Refs. [26 
27, S3]; for instance, in Ref. [28] the complete set of solutions for the Bogomol'nyi-Prasad- 



Sommerfield (BPS) states [25, 29 1 of the model was found. We notice that in the model 
described by W r (4>, x), the case r = lead us back to the standard 4 model. 

More recently, domain walls have been used in high energy physics as seeds for gen- 
eration of branes in scenarios where scalar fields couple with gravity in warped spacetime. 



In particular, in Ref. [15] we have used a model described by a single real scalar field to 



offer a scenario where the brane engenders internal structure. The model used in [15] was 
first introduced in Ref. |$3], and it depends on a real parameter, which generalizes the 
standard (ft 4 model, changing the way the scalar field self-interacts. It offers an interesting 
alternative to models which use two real scalar fields [10, Q as seeds for brane formation 



in warped spacetime in (4, 1) dimensions. The model is defined by [30 



P 



2p-l 



0(2p-i)/p 



P 



2p+l 



^(2p+l)/p 



(1.3) 



where p is an odd integer, p = 1,3, 5, ... The case p = 1 reproduces the standard 4> A model. 
The models defined by W p ((j)) support BPS solutions of the form 



4> p (x) =tanh p (x/p) 



(1.4) 



The case p = 1 reproduces the standard kink, but for p = 3,5, ... one gets new solutions, 



which very much resemble two- kink solutions [30]. In the above model, however, the scalar 
field self-interacts in an unconventional manner. For this reason, in the present work we 
return to the issue concerning thick branes that engender internal structure in models 
which support conventional self- interactions. 



As one knows [31], the presence of intricate self- interactions in general responds to 
thicken the brane pL fy|. Thus, in the present work we offer a way to trade the intricate 
self- interactions contained in the above model ([1.3]), to a more conventional model described 
by two real scalar fields, with the potential defined by Eqs. ( |1 . 1[ ) and (|1.2|), Within this 
context, it is interesting to see that the model defined by ( |1.2j ) is an extension of the 
standard i^> 4 model to the case of two real scalar fields, while the model defined by ( |l.3| ) 
extends the cfi model to the case of more intricate self- interactions. To implement the 
investigation, we split the work in the several sections that follow. In the next Sect. |2| we 
provide a complete investigation, which shows how both the Ising and Bloch walls translate 



- 2 - 



to stable branes in scenarios where the scalar fields couple to gravity in warped spacetime. 
In Sect. H we show how the second field required to implement the Bloch brane scenario 
contributes to engender internal structure to the brane solution. Finally, in Sect. ^ we 
present our comments and conclusions. 



2. Structure formation 

In this Section we first deal with domain walls in flat spacetime, which includes the basic 
issues one needs to investigate Bloch brane in curved spacetime. For clarity, we split the 
subject in two subsections, one dealing with Bloch walls, and the other with Bloch branes. 
We start the investigation with Bloch walls in flat spacetime in the subsection that follows. 

2.1 Bloch wall 

We first work in flat spacetime in (3, 1) dimensions. We consider the action 



d^x 



(2.1) 



for fj,,u = 0,1,2,3, with metric ds 2 = r) fll ,dx tJ 'dx u , where diag(ry Atl/ ) = (1,-1,-1,-1). We 
consider the potential 

V(<j>, X ) = ^(l-tf- r X 2 ) 2 + 2r W (2.2) 

which is defined by Eq. ( |1 . 1| ) , with W = W r ((j), x) given by Eq. (L2). This potential has 
minima at the points (±1,0) and (0, ±l/^/r) for r positive. 

The equations of motion for static solutions eft = 4>(x) and x = x( x ) that follow from 
the above model are solved by the first-order differential equations 

% = - 2r <t>x P' 4 ) 

We solve these equations using two distinct orbits connecting the minima (±1,0): first, 
we consider x = 0, which identifies a straight line segment joinning the minima (±1,0). In 
this case we get the one-field solution (f)(x) = tanh(x). Second, we consider the orbit 

<t> 2 + -^x 2 = i (2.5) 

for r G (0,1/2). This orbit identifies elliptical segments connecting the minima (±1,0). In 
this case we get the two-field solutions 

cp(x) = tanh(2rx) (2.6) 

X{x) = ±W~ - 2 sech(2rx) (2.7) 

We notice that the limit r — ► 1/2 changes the two-field solution to the one- field solution. 
The two-field solutions represent Bloch walls, while the one-field solution is the Ising wall. 
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We make this point clearer by recalling that the Ising wall is described by a straight line 
segment joinning the minima (±1,0) in the (4>,x) plane. For the Bloch walls, however, the 
orbits represent elyptical segments connecting the minima (±1,0), which may be used to 
describe chiral interfaces. The above model supports other solutions, as found for instance 
in Refs. |27|, [28 1. 

2.2 Bloch brane 

We now change to curved spacetime. We consider the action 



Sr. 



d xdy\/\g 



±R + \d a 4>d a ^ ± \d aX d a x - v{4>, x) 



where g = &et(g a b) and the metric is 

ds 2 = g ab dx a dx b = e 2A rj, lu dx fl dx 1 ' - dy 2 



(2.8) 



(2.9) 



where a,b = 0,1,2,3,4, and e 2A is the warp factor. We are working in (4,1) spacetime 
dimensions, using x 4 = y to identify the extra dimension. 

We suppose that the scalar fields and the warp factor only depend on the extra coor- 
dinate y. In this case the equations of motion reduce to the simpler form 



(/)" + AA!<j>' 
A" 



d(j> 



dx 



2 
3 



* 2 + x' 2 ) 



A 



1 



^ 2 ± x' 2 ) 



1 



V(<P,X) 



g K r ~ / 3 

where prime stands for derivative with respect to y. Our model is given by 



i 



dWr 



+ 



dWc 
dx 



-W 2 

3 c 



(2.10) 
(2.11) 
(2.12) 
(2.13) 

(2.14) 



where W c = 2W r , that is, it is twice the function we have used in flat spacetime; see, e.g., 
Eq. 0. 

The above potential leads to the following first-order differential equations, which also 
solve the equations of motion ||] 



1 dW c 

2 ~~d$ 

1 dW r 



and 



X 



A' 



2 8 X 



(2.15) 
(2.16) 

(2.17) 



We solve the first-order equations for <f> and x t° § e t the same results we have already 
obtained in flat spacetime, given by Eqs. (2.6) and (p. 7]), now changing x — ► y, to represent 
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the extra dimension. They are depicted in Fig. [1] for several values of the parameter r. 
We see that the width of these solutions depends on l/r, so it increases for decreasing r. 
We now turn attention to the behavior of the warp factor. We solve Eq. ( 2.17| ) to get 

A(y) = — (1 - 3r)tanh 2 (2ry) - 2 In cosh(2ry) (2.18) 

We plot the warp factor e 2A ^ for several values of the parameter r in Fig. [2]. It charac- 
terizes the presence of thick branes, in accordance with the behavior of the solutions (| 
and ( |2.7D for the scalar fields. 




Figure 2: Plots of the solutions of the warp factor exp[2 J 4(y)] for the values r — 0.05, 0.1 and 0.3. 

We now consider fluctuations of the metric and scalars. The perturbed metric is 

ds 2 = e 2A(y) + th ilv )dx ii dx v - dy 2 (2.19) 

and we consider <j> — > (f> + e(j) and \ ~~ * X + e X- We are using h^ u = h^ u (x, y), (ft = <f>{x, y), 
and x = xi x i V) to represent the respective perturbations; here x represents (x°, x , x 2 , x 3 ). 
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We consider variation of the action with respect to the scalar fields up to second order in 
e to obtain the equations for the scalar fluctuations. The procedure follows Ref. 32] and 

we get, for cf> 



e □0-4A0-0 +(f> hx d(j)dx = * 



and for % 



e U X 4,4 X X + X ^ 2 +0 - 2 xr? (_..!) 



where □ = d^d^. Also, we vary the action with respect to the metric up to second order 
in e to obtain 

-inv + e M Q«J + 2A'a„j / v - ^f'(d„d„h Xp - 8„9aV - MiV) 

+i*^"A'« W ^fc v ) + |e»V (i 8 -^ 1 + X^^) = (2.22) 

We choose as transverse and traceless ||, that is, we change h^ v — > h^ v = P^ v \ p h xp , 
where the projector is given by 

1 1 

PfivXp = -^nX^vp + K^plTvx) - -KpvTTXp (2.23) 



with 



= V,» - ^ (2-24) 



□ 

In this case the equation for fluctuation of the metric simplifies to 

h!' tlv + AA'h' lxu = e- 2A Uh llv (2.25) 

which describes linearized gravity, in the transverse and traceless sector. We notice that 
in this sector, gravity decouples from the matter field. 

We use another coordinate, z, which is defined by dz = e~ A<yV ^dy. Also, we set 

V(x, z) = e ik - x e~l A ^H^(z) (2.26) 

In this case the above equation ( |2.25|) becomes the Schrodinger-like equation 

d 2 H„ 



^ +U r (z)H^ u = k 2 H^ (2.27) 



dz 2 



where the potential is given by 



U r (z) = ^A"(z) + 9 -A' 2 (z) (2.28) 

One can show that the Hamiltonian in eq. ( 2.27] ) leads to no tachyonic excitations [15]. For 
the zero modes (k = 0) we have 

H^(z) = N^e 3 ^ 2 (2.29) 
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where N^u is a normalization factor. In Fig. [3] we plot the potentials U r (z) and the 
corresponding zero modes for some values of r. We see that U r (z) is a volcano-like potential, 
which may lead to the appearance of resonant states in its spectrum. For this reason, in 
Fig. [4] we depict some wave functions for massive modes in the two regions where k 2 < 
max[U r (z)] and k 2 > max[U r (z)]. The results show that the wave functions describe motion 
in the bulk, not bound to the brane, with no evidence for resonance when k 2 < max[[/ r (z)]. 
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Figure 3: Plots of the potentials Uq.o5(z), Uq.i(z), and (1/2)Uq.3(z) (right) and the corresponding 
zero modes (left). 




The above calculations show that the Bloch brane is consistent and cannot be discarded 
as an alternative to the brane scenario. More importantly, the Bloch brane scenario nicely 
leads to the possibility of brane with internal structure, as we investigate in the next 
section. 

3. Internal structure 

To see how the Bloch brane supports internal structure, let us now pay closer attention to 
the volcano potential depicted in Fig. [3] . There we see that it has a true minimum at y = 
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for 0.5 > r > r*, with r* ~ 0.17. However, when r reaches the interval r* > r > 0, the 
minimum at y = splits into two separate minima, indicating the appearence of internal 
structure. We make this point clearer with the help of our former work on brane structure 
1 15], from which we very naturally realise that the behavior of the volcano potential as r 
diminishes signals for the opening of an internal structure for the Bloch brane. 

To make the reasoning unambiguos, however, let us return to flat spacetime to focus 
on the issue in the simpler case firstly. We consider the system of Eqs. ( |2.3j ) and ( |2. -1| ) 
for the flat case to see that it decouples at r = 0, leading to (p(y) = tanh(y) in that case. 
However, in the limit r — > the two-field solutions ( |2.6| ) and ( |2.T| ) do not reduce to the 

are exact solutions 
0, so they are 



former solutions. The reason for this is that the solutions (2.6) and (£ 

for the system of first-order equations; they are not perturbative around r 
not expected to reproduce the solutions in the limit of vanishing r. Thus, it is exactly the 
nonperturbative behavior of the solutions that induce the appearence of nontrivial behavior 
for small values of r, for the specific nonlinearity of the system under consideration. 

To see how this works explicitly, we investigate the energy density of the nonpertur- 
bative two-field solutions in flat spacetime. We have 



1 /#(y)V 1 fd x (y)\ 



2 1 fdWrV 1 / dW r N 2 



+ 



(3.1) 



We use the solutions of the matter fields ( |2.6| ) and (|2 
we plot in Fig. [5] for several values of r. 



to obtain the energy density, which 
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Figure 5: The matter energy density in flat spacetime. The plots show £ a 05 (y),£ 01 (y), £ 0/25 (y) 
and£ a3 (y). 



We see that for 1/2 > r > 1/4, the energy density has the single-peak behavior around 
y = 0, as in the standard 4 model. However, for r* = 1/4 the plot shows the appearence of 
a plateau characterized by the maximum of dx/dy, signalling the splitting of the maximum 
into two new maxima, as depicted in Fig. [5]. Here we recognize that it is the x~ held which 
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responds for changing the behavior of the energy density for r in the interval 1/4 > r > 0. 
If we identify the center of the defect with the maximum of the energy density, we see 
that for r* > r > there appear two centers, identifying two interfaces which induce 
the appearence of internal structure. Moreover, the thickness of the defect increases for 
decreasing r, leading to the nice picture where as r decreases from r = 1/2, the defect gets 
thicker and thicker, until r reaches the interval (1/4,0), where two interfaces springs inside 
the defect. 

We now turn attention to curved spacetime, to the Bloch brane scenario investigated 
above. In this case, the matter energy density has the form 



1 /d^(y)\ 2 1 fd x (y)V 

2 V dy ) 2 V dy ) 



+ V c (cf>(y), X (y)) 



(3.2) 



We use the solutions of the matter fields and the warp factor (|2.6| ), and (|2.18| ) to 

obtain the matter energy density as a function of the extra dimension, which we plot in 
Fig. [6] for several values of r. 




Figure 6: The matter energy density in curved spacetime. 
E^{y),£^{y) & nd{l/2)E^(y). 



The plots show £° m (y), 



The investigation shows that for 0.5 > r > 0.17 we have the single-peak behavior in 
the very inside of the brane. For r* ~ 0.17 we also find a plateau, as we did in the flat 
case. Here, however, it appear at the lower value r* < r* = 0.25, due to the coupling to 
gravity. But in the region 0.17 > r > the behavior for the energy density shows the 
appearance of two peaks, as we illustrate in Fig. [6]. Thus, the coupling to gravity does not 
destroy the behavior which induces the presence of internal structure for the Bloch brane. 
Moreover, the same effect of increasing the brane thickness as r goes to zero is observed. 
Furthermore, we emphasize that for the same r, the thickness of the defect diminishes in 
the curved case, in comparison with the flat case. This is compatible with the reasoning 
that gravity induces an attraction between the interfaces, as we have already identified in 
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Ref. [15] — compare Figs. [5] and [6]. Since the value of r* which induces the appearance 
of internal structure is lower in curved spacetime, we can infer that gravity tends to supress 
the presence of internal structure in braneword scenarios like the one here investigated. 



4. Further Comments 

In this work we have investigated the presence of Bloch brane in braneworld scenarios 
which appear in models involving real scalar fields coupled to gravity in (4, 1) dimensions 
in warped spacetime. Our investigations are based on several recent works on branes 
0, |, § |, H|, 0, |, |, and it follows the lines of our former work on brane structure |15| . 



Former investigations on braneworld models described by scalar fields coupled to grav- 
ity in warped geometry such as the ones in Refs. g, |5| show that the presence of scalar 



field with intricate self-interactions directly contributes to thicken the brane. Furthermore, 
in [15] one shows that self- interactions may also lead to brane engendering internal struc- 



ture, in a way similar to that found in Ref. [10]. This fact has given the main motivation 



for the present work: to show how the unconventional self- interactions introduced in [15] 
can be traded to more conventional interactions, under the expences of adding more fields 
to the matter field contents of the model. The Bloch brane solution here introduced is 
stable, and it correctly responds for the presence of nontrivial structure inside the brane. 
The addition of the second field allows to trade intricate self-interactions present in the 
single-field scenario to the conventional interactions that appear in the two-field model. 
A natural extension of the two-field model, leading to three fields could be 

W(^x,p) = ^-\^-r<i ) {x 2 +p 1 ) (4.1) 

which was already considered before in other contexts J33| , |34] ]. In this case the first-order 
equations for the matter fields are solved by 

<p(y) = tanh(2ry) (4.2) 

and 



\ ! .r ) - \/--2 cos(0)sech(2r?/) (4.3) 



p(x) = J - - 2sin(6>)sech(2ry) (4.4) 

with 1/2 > r > 0, and 2ir > 6 > 0. The new parameter 9 drives the elliptical orbits to form 
an ellipsoid. This extension changes the former scenario with the inclusion of a new degree 
of freedom, which may be used to make the second field complex, allowing for the presence 
of the global symmetry U(l), which can be gauged to lead to more realistic models. Other 
extensions to three-field models could be done. An interesting line of investigation could 



follow Ref. [35 1, which investigates how walls and global and local strings can end on walls. 
Such investigations would require the presence of two or more extra dimensions. This and 
other related issues are presently under consideration. 
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